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In the Ginzburg-Landau theory of superfluid "^He, the free energy is expressed as an expansion 
of invariants of a complex order parameter. Strong coupling effects, which increas e w ith increasing 
pressure, are embodied in the set of coefficients of these order parameter invariant^^. Experiments 
can be used to determine four independent combinations of the coefficients of the five fourth order 
invariants. This leaves the phenomenological description of the thermodynamics near Tc incomplete. 
Theoretical understanding of these coefficients is also quite limited. We analyze our measurements 
of the magnetic susceptibility and the NMR frequency shift in the B-phase which refine the four 
experimental inputs to the phenomenological theory. We propose a model based on existing ex- 
periments, combined with calculations by Sauls and Seren^ of the pressure dependence of these 
coefficients, in order to determine all five fourth order terms. This model leads us to a better un- 
derstanding of the thermodynamics of superfluid ^He in its various states. We discuss the surface 
tension of bulk superfluid ''ife and predictions for novel states of the superfluid such as those that 
are stabilized by elastic scattering of quasiparticles from a highly porous silica aerogel. 



I. INTRODUCTION 

The Ginzburg-Landau (GL) formulation gives a phe- 
nomenological representation of the free energy of su- 
perfluid " ^He as an expansion in terms of the order 
parametei!32l3|. The expansion coefficients specify the 
stability of various p-wave states and their thermody- 
namics near Tc- These coefficients are well-defined the- 
oretically for the weak coupling case. However, '^He is 
not a weak coupling superffuid as is clear from its phase 
diagram where there is a region of ^-phase at high pres- 
sures. This is in contrast to the weak coupling limit for 
which the B-phase is always stable. The strong cou- 
pling correction to the pair interaction is responsible for 
the A-phase, an effect of spin and density ffuctuations 
proportional to Tc/T^. Calculations'^ cannot account 
quantitatively for the strong coupling corrections and so 
the coefficients must be determined empirically. Five of 
these parameters are coefficients of the fourth order in- 
variants of the order parameter in the GL free energy. 
These are called the /3-parameters, /3i's where i = 1, 
5. Unfortunately, there are not enough independent sets 
of experiments to determine all the parameters and so the 
phenomenological description of superffuid ^He is under 
determined. This hampers our ability to predict stabil- 
ity for novel superffuid p-wave states, such as those that 
might be favored by elastic scattering from high porosity 
silica aerogel. 

In this paper, we present four combinations of (3iS 
which we determine from measurements and we describe 
a model which resolves the ambiguity in identifying all 
ffve of them independently. The coefficient of a ffeld de- 
pendent term in GL theory, gz , plays an important role in 
determining more accurate combinations of the [3^ than 
have been previously reported. Our NMR measurements 
of the susceptibility^] show that Qz is close to its weak 
coupling value at all pressures. This allows us to in- 
terpret our high resolution measurements of the NMR 



frequency shift in the B-phas^^ and to obtain accurate 
/3-parameter combinations. 

Our model for determining the five/3i 's is motivated by 
the calculations of Sauls and Serene^. We note that the 
calculations, although only qualitatively consistent with 
the existing experiments, nonetheless can accurately ac- 
count for their pressure dependence. Furthermore, we 
note that the calculations indicate that one of the /3- 
parameters, is close to its weak coupling value at 
all pressures. Motivated by these observations and the 
fact that the experimentally known combinations of the 
/3i's approach their weak coupling values at zero pres- 
sure to within 5%, we make the following two assump- 
tions: First, the /3-parameters are, on average, close to 
their weak coupling values at zero pressure and we use 
this criterion to select (3i at zero pressure. Second, we 
take their pressure dependences from the theory which 
seems to accurately represent this aspect of the known 
/3-parameter combinations. These assumptions are suf- 
ficient to constitute a model to determine the full suite 
of /3-parameters. With this information we can calculate 
the surface tension between A- and B-phases in bulk su- 
perfiuid "^He and compare with experiment. We can also 
calculate the stability of the axi-planar state in bulk su- 
perfluid ^He as a function of pressure and we can evaluate 
predictions for anisotropic p-wave states that are robust 
in the presence of elastic scattering from silica aerogel. 



II. GL THEORY FOR ^HE 

A phenomenological macroscopic description of phase 
transitions is given by the GL theory, in which the free 
energy is expressed as an expansion of the order parame- 
ter. In the case of superfluid ^He, the order parameter'^, 
A, is a complex 3x3 matrix and the free energy of the 
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system can be expressed as, 

F = ~aTY{AA^) + g,H^{AA^)^,H, + f3,\TY{AA'^)\'' 
+/32[Tr(AAt)]2 + /33Tr(A^^(AA^)*) 
+/34Tr((AAt)2) + (3r,Tr{AA^AA^)*). (1) 

Here the dipole energy term is neglected. The magnetic 
field components are and A^ and A'^ are the Her- 
mitian conjugate and transpose of A. The structure of 
the order parameter admits five fourth order invariants 
each of which has a corresponding coefficient, /3i. At the 
second order thermodynamic transition to superfiuidity, 
Tc, all p-wave superfluid states are equally probable, but 
their stability below depends on the Pi. In the weak- 
coupling limit the free energy coefficients are. 
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where the normal density of states at the Fermi energy 
is N{fi), the gyromagnetic ratio for '^He is 70, fcs is the 
Boltzmann constant, Fq is a Fermi liquid parameter de- 
termined from the magnetization measurement!^ and C,{x) 
is the Riemann zeta function. However, '^He is not a 
weak coupling superfiuid and strong coupling effects in- 
crease with pressure. The strong coupling corrections 
for a are negligible but they have a significant effect on 
the PiS, and might also contribute to gz- Calculations 
of strong coupling corrections have been performed for 
model potentials^ those of Sauls and Serene^J being the 
most complete and the ones we will refer to in this work. 



III. EXPERIMENTS 

There are seven free energy coefficients which must 
be determined from experiment. The difficulty lies in 
the fact that there is insufficient experimental input to 
constrain this phenomenological description of superfiuid 
•^He. Among the seven coefficients, a and gz are de- 
termined without ambiguity. The measurements of the 
specific heat in the normal state, Cjv and the transition 
temperatur^i^ Tc, give us a. The slope of the "^He-i? 
magnetizatiorP extrapolated to dMs / dT\T^, and the 
specific heat jumjJI^ of '^He-i?, ACb/Cn, are required 
for gz, for which we have new results presented in this 
section. For the remaining five /3i's, there are only four 
independent sets of experiments so that only four combi- 
nations of Pi's can be found in the form of sums. These 
are P345, P12, P245, and P^, where we use the Mermin- 
Stare convention, Pij — Pi + Pj. 

First, we will describe the relevant experiments and 
the logic for determining these combinations of the Pi's. 



A) . P345 requires measurements of the ^He-B trans- 
verse NMR (7-shiftP^, g, which must be combined with 
the slope of th e B-phase longitudinal NMR resonance 
frequencjdSEIl^ — t), in the limit approaching 

Tc as well as with measurements of ACb / CjJ^ where 
t ~ T/Tc. In order to have the value of the g-shift at Tc 
it is helpful to observe that the i3-phase susceptibility and 
the g-shift are linearly related, facilitating an extrapola- 
tion to Tc. The B-phase heat capacity jump is measured 
only below the polycritical point (PCP). However, mea- 
surement of the specific heat in the A-phase along with 
measurements of the latent heat at the A- to B-transition 
allows a thermodynamic calculatiorP^of the specific heat 
jump for the B-phase at pressures above the PCP. Con- 
sequently, ACb/Cn is experimentally determined at all 
pressures. 

B) . From the specific heat jump, ACb/Cn and the 
values for P34C, obtained above we can directly determine 

Pl2- 

C) . From the specific heat jump, ACa/Cn we can di- 
rectly determine P245 , but only for pressures greater than 
the PCP where this jump can be measured. Below the 
PCP P245 is found from the quadratic magnetic field sup- 
pression of the first order '^He A- to B-transitiorP^, g{P), 
along with the values of P12 and P345 that have been 
obtained above in A) and B). 

D) . Finally, we can fix P^^ uniquely by the asymmetry 
ratio, r, of the linear field dependent splitting of the Ai 
to A2 transitions^'* in high magnetic field combined with 

P245- 

In summary, four independent combinations of ex- 
periments gives us four constraints on the Pi's, which 
is insufficient to identify all five of them. In princi- 
ple, measurement of the surface tension at the '^He A- 
B interface could provide us with a fifth independent 
combinatiorfiSESE^ of pi's. However, the surface tension 
vanishes near Tc due to the degeneracy of the free energy 
at Tc of A- and B-phases. For this reason it is not pos- 
sible to obtain sufficiently high resolution measurements 
of the surface tension to provide useful characterization 
of strong coupling effects in the Ginzburg-Landau limit. 
In the following, we will discuss in more detail the exper- 
imental determination of strong coupling and its effects 
on the Pi's. 

The coefficient for the field coupling term, gz, is de- 
termined by measuring the slope of the magnetization of 
•^He-i? in the limit approaching Tc, 
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dm 
dt 

' dm 



AC^" 
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where m = Mb/M^ and Mjv is the normal state mag- 
netization. The superscript wc, which we use here and 
in the following, indicates the weak coupling limit. 

Magnetization measurements of superfluid '^He 
have been of great interest since its discovery. 
Two differe nt techniqu es - NMR based dynamic 
measurementiPiSliniinili] and SQUID based static 
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TABLE I: Ginzburg-Landau /3-parameters and the experimental quantities from which they are derived. The NMR j?-ph ase 
g-shift is a fit to data from Kycia^ given by HaarcP'. The NMR _B-phase longitudinal resonance was measured by RancP^^H for 
which a smoothed fit is given by Haard^. The coefficient of quadratic magnetic field suppression of the _B-phase was measured 
by Tang et al^^. The B-phase heat capacity jump was taken from GreywalP^ and the asymmetry ratio of the linear field 
dependent splitting of the Ai to A2 transitions was reported by Israelson et aZ.'^. Extension of the measuerments of the A- 
phase heat capacity jump to pressures lower than the PGP requires a calculation based on the measured quadratic suppression 
of the A- to B-transition as described in the text. 



measurement^SlIIIMIlIl _ have been performed over 
the past thirty years. Historically, there has b een 
a discrepancy between these two technique^^MIl the 
origin of which has not been estabhshed. Nonetheless, 
more recent experiment d^ l ^^ l bring the results closer 
together. Haard measured the magnetization using high 
resolution NMBP. A careful analysis of this and other 
measurement42S122l reveals that the discrepancy appears 
to be negligible near the transition temperature Tc- 
Using Eq. 6, HaarcP determined gz from NMR and 
found the results presented in Fig. 1, where gz is close 
to its weak coupling value, i.e. gz — 1- From Haard's 
measurement, the deviation from weak coupling appears 
to grow slightly with pressure. F rom a nalysij^ of the 
more accurate work of Scholz et al^^^^ it appears that 



gz is pressure independent. The difference between the 
data sets is likely due to the wider range of extrapolation 
in the _B-phase toward Tc that is required to determine 
gz at elevated fields in the case for Haard's measurement. 
Hahn et aZ.l^came to the same conclusion, gz — 1, based 
on their SQUID measurements, and so we will take gz 
to have its weak coupling value at all pressures. Having 
established Pz, /3345 can be calculated from the NMR 
ff-shif t'^'^''''^ of the transverse NMR frequency in ^He-i? 
which has the following relationship!^ with cjz and : 

/3345 ^ /?3T5 ( Cn \ ''m ( h y 1 

gz [l + F^y^ \ ACb) l-t \2TTkBTj g ^ ' 
In earlier report^ of the g-shift, the analysis to ob- 
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FIG. 1: (jz obtained from magnetization measurements by 
NMR. Closed circles a re the measurements by HaarcP^and 
open circles by ScholJ^^'^. The resuhs from both mea- 
surements are consistent and give approximately unity 
for 



tain /3345 estimated gz incorrectly. The values in Table 
1 for the g-shift and the B-phase longitudinal resonance 
frequency are smoothed values^ from a large number of 
experiment^, significantly more than what was originally 
reported by Kycia et alW. GreywalP^ has measured the 
specific heat of "^He-A and B. The specific heat jump at 
Tc, for these two phases, is related to Pa and (3b through: 
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where a' = da/dT. At pressures less than the PCP, 
the magnetic suppressiorP^I, g{P), of the AB transition 
temperature, T^b, is used to obtain (^245 through: 



.^l + (fe//?^-l)(l + T3^) + l 
n/3) = TTll^ — • (10) 



Pb/Pa - 1 



Here g{P) is defined by. 



Tc -^^^^ \Bo 




(11) 



where B is the applied magnetic field and Bq = 
N{0)/6gz- Finally, /Js can be determined by measuring 
the asymmetry ratid^ilof the A1-A2 splitting, r. 
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The four experimentally determined /^-coefficient com- 
binations, along with the measurements used to obtain 
them, are tabulated from to 34 bar in Table I. 



IV. MODEL FOR DETERMINING (3'S 

As stated earlier, we impose two assumptions to elimi- 
nate ambiguity associated with sorting out all five Pi's 
from the four known combinations of Pi's determined 
from the experiments described in the previous section. 
The assumptions are: 1) the pressure dependence of Pi 
calculated by Sauls and Serene^ is valid. 2) At zero pres- 
sure, all five Pi's approach their weak coupling values, on 
the average. The consequences of these assumptions will 
be discussed in the following subsections. 



A. Comparison with the Calculation 

Sauls and Seren^ developed a potential scattering 
model to find the strong coupling corrections to the 
/3-coefficients in the pressure range of 12 to 34.4 bar. 
Since we do not have five experimentally determined P- 
coefficients with which to directly compare to the theory, 
we construct from the calculation those four combina- 
tions of /3-coefficients, /?345, P12, P245, and P5 that are ex- 
perimentally accessible and compare these with the mea- 
surements in Fig. 2. First, we note that the experimental 
results suggest that superffuid ^He is predominantly weak 
coupling at zero pressure. Secondly, the pressure depen- 
dence of each combination shows remarkable agreement 
between experiment and theory for P > 12 bar, the range 
where the calculations were performed. It is also appar- 
ent that the calculation of the absolute values of the Pi's 
is less reliable than their pressure dependence. Finally we 
note that, in the calculation, the smallest strong coupling 
correction among the Pi 's is for Pi . Guided by this infor- 
mation, we will assume that the pressure dependence of 
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FIG. 2: Comparison of four known /3-combinations from 
the experiments (dashed lines) and Sauls and Serene's 
calculatiori-^ (solid lines). The pressure dependences are 
in good agreement but the absolute values are not as 
close. 
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Pi (P) can be taken from the Sauls and Serene calculation 
and then we need only determine /?i(0). 



V. APPLICATIONS 



Surface Tension at the A-B Interface 



B. 



Zero Pressure Values of the /3i's 



The pressure dependence of the (3i's is insufficient to 
resolve the ambiguity associated with the f3 coefficient 
combinations. Five independent values of /3i's at a given 
pressure are required along with the pressure dependence 
for /3i. The calculations indicate that strong coupling 
corrections are smallest for Pi{P) and from experiment 
we see that the measurable combinations deviate from 
their weak coupling values by less than 5% at zero pres- 
sure. On this basis one possibility would be to simply 
choose f3i{0)/f3Y^ = 1, i.e. to be weak coupling. Another 
possibility, the more democratic one, is to choose /?i(0) 
as a variational parameter and minimize the mean square 
deviations of all /9-parameters from their weak coupling 
values at zero pressure subject to the constraints imposed 
by the four different combinations that have been deter- 
mined experimentally. For the latter method we find 
Pi{0)/PY'^ — 0.97 which is essentially equivalent to the 
first choice. In the following we make the latter choice. 
We show this process explicitly in Fig. 3 where we calcu- 
late all of the /3i{0)'s as a function of /3i(0) subject to the 
four experimental constraints. It is clear that for /9i(0) 
near its weak coupling value, as emphasized by the cir- 
cled region, all the others approach their weak coupling 
values at zero pressure as well. With this choice for /3i(0) 
and the pressure dependence of (3i taken from Sauls and 
SerentP, is now uniquely defined and all the other 

Pi's can be determined. These Pi's are tabulated in the 
first five columns of Table II. 
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With all the Pi's now determined, one can calculate 
the surface tension between the A- and _B-phases of su- 
perftuid ^He. According to Thuneber^Sl, the surface free 
energy of the A-B interface is expressed as 
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where ^(T) is the temperature dependent coherence 
length of ■^He, 
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1.89 - 1.98V^ - 0.31k for k < 1/V2, 
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Here a and c are defined as a = 2/3i ~V P-i — /345 and 
c = -(2/3i + /3345). The pf^'s are any set of A's that 
satisfy the condition for the surface energy to vanish, 
2/?i 4-/33 = 0,/345 = 0. Weak coupling values of Pi's are 
a subset of the P^^'s, but the P^l'^'s need not be limited 
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FIG. 3: Zero pressure values of the /3i's parameterized 
by /3i(0). The numbers on each line correspond to the 
subscript i oi Pi. The /3i(0)'s clearly converge around 
/9i(0)//3™'^ = 1, which is an indication that the /3i's tend 
toward their weak-coupling values at low pressure. 



FIG. 4: Osheroff and Cross's measurements of sur- 
face tensiorP^ in comparison to the calculation by 
Thuneberg^ for various /3i choices and Ijpf^ ~ /9245 + 
(/3i2 + /3345 / 2)~^. The comparison of the two with the 
choice of Pi/ P\'^ ~ 1 is consistent, but the measurement 
lacks the resolution to conclusively determine the value 
of Pi. 
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-0.20 


0.51 


0.36 


0.48 


-0.50 


7 


-0.98 


1.98 


1.70 


1.44 


-1.82 


-0.28 


0.72 


0.53 


0.66 


-0.70 


8 


-0.98 


1.98 


1.70 


1.42 


-1.83 


-0.35 


0.87 


0.66 


0.78 


-0.84 


9 


-0.99 


1.98 


1.69 


1.41 


-1.84 


-0.41 


0.99 


0.75 


0.87 


-0.95 


10 


-0.99 


1.97 


1.69 


1.40 


-1.85 


-0.45 


1.08 


0.83 


0.94 


-1.05 


11 


-0.99 


1.97 


1.70 


1.39 


-1.86 


-0.49 


1.15 


0.90 


0.99 


-1.12 


12 


-0.99 


1.96 


1.69 


1.39 


-1.87 


-0.52 


1.21 


0.96 


1.03 


-1.17 


13 


-0.99 


1.95 


1.69 


1.39 


-1.88 


-0.55 


1.26 


1.01 


1.06 


-1.22 


14 


-1.00 


1.95 


1.70 


1.38 


-1.89 


-0.58 


1.30 


1.05 


1.09 


-1.27 


15 


-1.00 


1.95 


1.72 


1.35 


-1.89 


-0.60 


1.34 


1.10 


1.10 


-1.32 


16 


-1.00 


1.95 


1.73 


1.34 


-1.90 


-0.62 


1.38 


1.13 


1.12 


-1.35 


17 


-1.00 


1.94 


1.72 


1.33 


-1.90 


-0.64 


1.40 


1.16 


1.13 


-1.38 


18 


-1.00 


1.94 


1.73 


1.32 


-1.91 


-0.66 


1.43 


1.19 


1.14 


-1.41 


19 


-1.00 


1.93 


1.72 


1.33 


-1.92 


-0.67 


1.45 


1.21 


1.16 


-1.43 


20 


-1.01 


1.94 


1.74 


1.31 


-1.93 


-0.69 


1.48 


1.24 


1.16 


-1.47 


21 


-1.01 


1.94 


1.74 


1.29 


-1.93 


-0.71 


1.50 


1.26 


1.16 


-1.49 


22 


-1.01 


1.93 


1.74 


1.29 


-1.94 


-0.72 


1.51 


1.28 


1.17 


-1.51 


23 


-1.01 


1.93 


1.74 


1.29 


-1.95 


-0.73 


1.53 


1.30 


1.18 


-1.53 


24 


-1.01 


1.93 


1.74 


1.28 


-1.96 


-0.74 


1.54 


1.32 


1.18 


-1.54 


25 


-1.01 


1.93 


1.74 


1.28 


-1.97 


-0.75 


1.56 


1.33 


1.18 


-1.58 


26 


-1.02 


1.93 


1.73 


1.27 


-1.97 


-0.76 


1.57 


1.34 


1.18 


-1.60 


27 


-1.02 


1.93 


1.74 


1.26 


-1.98 


-0.77 


1.58 


1.36 


1.18 


-1.61 


28 


-1.02 


1.93 


1.73 


1.26 


-1.99 


-0.78 


1.60 


1.37 


1.19 


-1.62 


29 


-1.02 


1.93 


1.73 


1.26 


-2.00 


-0.78 


1.61 


1.38 


1.19 


-1.63 


30 


-1.02 


1.93 


1.72 


1.26 


-2.01 


-0.79 


1.62 


1.38 


1.19 


-1.67 


31 


-1.03 


1.93 


1.73 


1.25 


-2.02 


-0.80 


1.62 


1.40 


1.19 


-1.68 


32 


-1.03 


1.93 


1.73 


1.25 


-2.02 


-0.81 


1.63 


1.40 


1.19 


-1.68 


33 


-1.03 


1.93 


1.73 


1.25 


-2.03 


-0.81 


1.63 


1.41 


1.19 


-1.69 


34 


-1.03 


1.93 


1.73 


1.25 


-2.03 


-0.82 


1.64 


1.42 


1.20 


-1.70 



TABLE II: (3i's for bulk superfluid "^He, left side, and superfluid '^He in 98% porosity aerogel in the IISM with A= 150 nm and 
^a=40 nm, right side. 



to their weak coupling values. Thuneberg suggested two 
different values for 2//?^°^ = /S^/g + (/?i2 + /3345 /2)-^ 
and f]^'' = Pi2 + Psisf^, but kept the relative magnitude 
of the five p'f^'^s the same as for the weak-coupling case 
in his original work^. We examined these two choices of 

From a number of different choices for the /3i including 
the values of the /3i chosen from Table I, the surface ten- 
sion at the melting curve is calculated. The calculation 
with our choice of /3i and the measurements of Osheroff 
and Crosd^ are in good agreement . An example of the 

calcuation with 2/l3^^ ~ /3245 + (/3i2 + /^345/2)^"^ is shown 
in Fig. 4. The calculation, however, has a number of 
limitations. One is that the calculation depends on the 
choice of which is not uniquely defined. The other 
is that the experimental results do not have high enough 
resolution to determine the /3i's independently. 



B. How Stable Is the Axial State? 

A number of experiments performed to investigate the 
order parameter of ^He-A phase have confirmed that the 
A-phase is, in fact, the axial state. This confirmation 
could be further strengthened by studying the thermo- 
dynamic stability of the axial state over other possible 
equal spin pairing states, such as an axi-planar state; 
some concern has been raised in the past that an axial 
state and an axi-planar state may not be easily distin- 
guishable due to the ir continuously related order param- 
eter structure^^Ml However, a certain combination of 
/3-coefficients, namely P45 can be used to check the rela- 
tive thermodynamic stability between the two states. If 
P45 is negative, the ^-phase is the axial state and if P45 is 
positive, the A-phase would be the axi-planar state. By 
imposing f]^^ = as the fifth constraint in addition to the 
four known combinations of the f3i's, a unique set of Pi's 
is obtained which can be used to plot a phase diagram 
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FIG. 5: Phase diagram for axial and axi-planar states 
with /3i as a parameter. The choice of /3i with our model 
(dashed line) places the ^-phase in the region of the axial 
state. 



for axial and axi-planar states with f3i as a parameter. 
This phase diagram is shown in Fig. 5. We compare 
our choice of /3i from Table I with the phase diagram 
and this value lies well within the axial state regime at 
high pressure as has been commonly believed and which 
a number of experiments independently confirm^^ 1233^ 
However, it should be noted that our choice of (3i indi- 
cates that there is a near degeneracy of the axial and 
axi-planar states at zero pressure and this might be in- 
teresting to investigate further. 



C. The Robust Phase in Aerogel 

Full determination of all five f3i 's has important impli- 
cations for superfluid ■^He in aerogel. Within the con- 
text of scattering models we can calculate the appropri- 
ate modifications to the /3i's and explore the predicted 
stability of various superfluid states. There are multiple 
different scattering models for ^He in aerogel, e^ the ho- 
mogeneous isotropic scattering model (HISM)B2 and in- 
homogeneous isotropic scattering models (IISMp ^ l '^^ l ^^ l 
We use the IISM of Sauls and SharmaP^, a modification 
of the HISM of Thuneberg et al^. The Pi's in the IISM 
are modified through: 



f 0i\ 

Pi 

PI 



/ \ 



V-1/ 



PS 



V-2/ 
30(7rfcBTe)2 ^ (2n-l + 



( A/3f \ 



,(17) 



1 



(18) 



9(7rfcBT,)2 



- (2n 



X 



(19) 



where x = x/{l + Cl/x), (a = x = hvp/^TikBTX, 

is the strand-strand correlation length, A is the transport 
mean free path for "^He quasiparticles, and is the s- 
wave scattering phase shift. 

With the five PiS for bulk superfluid given in Table II, 
we calculated the effects on the Pi 's of scattering from the 
aerogel strands. We distinguish these coefficients from 
bulk ^He with a superscript, /3f . We assumed unitary 
scattering, = 7''/2 with A = 150 nm and = 40 nm. 
These parameters are typical of 98% porosity aeroge 
The effects of scattering in the weak coupling approxi- 
mation are included in both /3q and b. In addition, the 
Pf's will have a strong coupling component that will be 
modifled by elastic scattering. We accomodate this by 
rescaling the A^f^'s with a factor Tca/Tc, since strong 
coupling effects^l are linear in T^/Tp. The results of the 
calculation, Pf /Pq, are tabulated in the last flve columns 
of Table II. For this choice of aerogel parameters the su- 
perfluid state is not stable below a pressure of 5 bar as 
reported by Matsumoto et al.^^, and hence the table is 
blank below this pressure. 

A direct consequence of the modification of Pf's ac- 
cording to the scattering model is the enhancement of rel- 
ative stability of the _B-phase with respect to the A-phase 
for ^Hc in aerogel. For either the HISM or IISM, the 
isotropic state (i?-phase) is found to be stable over the 
entire pressure range. However, superfluid '^He in aero- 
gel has a metastable A-like phase that has been clearly 
observecPSliMl various samples on cooling below Tc- 
Although the exact nature of this phase is still in ques- 
tion, it is known that the metastable phase is an equal- 
spin-pairing state''^, similar to the bulk A-phase, hence 
it is referred to as an A-likc phase. However, lack of un- 
derstanding of the orbital part of the order parameter 
makes the identity of the state less clear. Furthermore, 
the question of stability of any equal-spin-pairing state 
with repect to the aerogel B-phase relies on an under- 
standing of the appropriate /3-parameters for which we 
have no direct independent information. VolovilP^l has 
argued that the axial state in the presence of quenched 
anisotropic disorder cannot exist as a spatially homo- 
geneous superfluid owing to arguments from Imry and 
MaP^. If the metastable phase is in fact the axial state, 
the order parameter would not have long range orien- 
tational order, a state which Volovik has called a su- 
perfluid glass. With a different approach, Fomin^l has 
argued that there are other p-wave pairing states which 
are also equal-spin-pairing but do not suffer from the 
same difficulty, and that these might be candidates for 
the metastable aerogel phase. Such phases would be 
robust in the presence of anisotropic scattering, mean- 
ing that A^iA*^^ + Af^jA*^^ cc (5,^ where S^j is the Kro- 
necker delta.HS NMR experiments have been performed 
on "^He in 97.5% aerogel which support the view that 
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FIG. 6: The asymmetry ratio, r, for the Ai — A2 spht- 
ting was calculated for the axial state (solid red line) and 
the robust phase (dashed line) in aerogel, where we used 
the IISM of Sauls and Sharma^^ with a transport mean 
free path for '^He quasiparticles A = 150 nm and strand- 
strand correlation length = 40 nm, which match well 
to phase diagram measurements on the same sample by 
Gervais et aZ.™ The measurements of H.C. Choi et al^^ 
(closed circles) are more consistent with the A-like phase 
of aerogel ^He being the axial state than the robust state. 



the metastable A-like phase is in fact a robust stat^^, 
but other measurement^^^llllll appear to be inconsistent 
with this interpretation. The free energy for the robust 



state^ can be expressed as, 



where. 



(/3i3 + 9/32 + 5/345 )/9. 



(20) 



(21) 



Thermodynamic properties of the robust state have 
not been predicted because it involves all five piS be- 
yond the four combinations known to us so far. However, 
the determination of j3iS from our model allows us to in- 
vestigate the properties of the robust state. First, we 
calculate the asymmetry ratio of the A1-A2 splitting in 
aerogel. For the yl-phase this ratio is expressed in terms 
of the Pi's given by Eq. 12. In the case of the robust 
state the ratio rn is given bj^^^^. 



rR 



15 



/3i3 + 9/32 + 5/345 ■ 



(22) 



With the values of the Pi's from Table II, the asymmetry 
ratio rn is found to be ~ 0.2, considerably smaller than 
what has been found experimentally^^, > 1.0. These 
results are compared in Fig. 6. 

Second, we calculate the relative stability of the robust 
state with respect to the B-phase over the pressure range 
from zero to 34 bar with Pi for bulk ^He as a parame- 
ter subject to the constraints of the four experimentally 



ca 



-2 - 





1 1 1 1 1 

from this worl< 


'3 


Isotropic 


LL 




irma 





z 


1 Robust 




1 1 1 1 1 



0-3 



-1 - 



10 15 20 25 

Pressure (bar) 



30 



35 



FIG. 7: Phase diagram for the isotropic phase and the 
robust phase for '^He in 98 % porosity aerogel with /3i 
(bulk) as a parameter. Our choice of fi\ (dashed line) 
makes the isotropic phase more stable than the robust 
phase. 



known combinations of the P's given in Table I. These re- 
sults are shown in Fig. 7. For the robust state of '^He in 
aerogel to be stable. Pi would have to be significantly dif- 
ferent from the value derived from our model, assuming 
the form of the free energy in Eq. 1. 



VI. CONCLUSIONS 

We have investigated the experimental basis for deter- 
mining strong coupling as a function of pressure for su- 
perfiuid '^He based on analysis of our NMR data. Given 
the limitation that we have only four experimentally 
identifiable /3-coefficient combinations, we developed a 
phenomenological model with two assumptions: 1) that 
superfluid ^He is predominantly weak coupling at low 
pressure and 2) that the pressure dependence of Pi can 
be taken from Sauls and Serene's calculation. This model 
provides us with all five /3-coefficients. Using this model 
we calculated the surface free energy at the A-B interface 
and compared with experiment. Although the measure- 
ment does not have high enough resolution to validate 
our model, it is not in disagreement. The model is also 
consistent with the general consensus that the so-called 
A-phase is the axial phase rather than the axi-planar 
phase. We used our values of the Pi's to calculate the 
corresponding strong coupling effects for superfluid '^He 
in aerogel. We find that the B-phase is stable at all pres- 
sures. We compared the relative stability of the robust 
state proposed by Fomin with that of the S-phase. The 
robust state is unstable relative to either the isotropic 
_B-likc phase or the axial state. Furthermore, the asym- 
metry ratio, r/j, of the Ai-A^ splitting for superfiuid 
•^He in aerogel was calculated for the robust state and 
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it was found to be significantly smaller than the experi- 
mental values. Our interpretation is that the A-like aero- 
gel phase is not a robust state based on the free energy 
expansion given in Eq. 1. 
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